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Abstract 

We study the consequences of including parity preserving matter 
for the effective dual theory corresponding to compact QED^; in par- 
ticular we focus on the effect of that contribution on the confinement- 
deconfinement properties of the system. To that end, we compare two 
recent proposals when massless fermions are included, both based on 
an effective anomalous dual model, but having global and local Z2 
symmetries, respectively. 

We present a detailed analysis to show that while for large mass 
fermions the global Z2 symmetry is preferred, in the massless fermion 
case the local Z2 scenario turns out to be the proper one. 

We present a detailed discussion about how the inclusion of mass- 
less fermions in compact QED3 leads to deconfinement, and discuss 
the stability of the deconfined phase by introducing a description based 
on an instanton dipole liquid picture. 
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1 Introduction 



Three-dimensional non Abelian gauge theories have many properties that 
make them very interesting from the point of view of the realization and de- 
scription of the confinement mechanism. Among those properties a remark- 
able one is that, in the low-energy regime, there exist effective descriptions 
which exhibit confinement in a very simple way, already in a semiclassical 
approach. As explained by 't Hooft [I], for the case of an SU(N) gauge field 
theory, the natural degrees of freedom of the effective field theory turn out 
to be magnetic vortex fields, with the confinement mechanism being tanta- 
mount to the confinement of topological solitons for those vortex fields. Using 
this approach, the phase structure of the SU(N) Georgi-Glashow model in 
2 + 1 dimensions can be obtained by showing that the magnetic vortices 
condense, as in this situation the heavy electric charges are confined. Indeed, 
Cy, the effective Euclidean Lagrangian for the dual theory is given by 

C v = d^Vd.V + X(VV - /i 2 ) M + c(Y N + V N ) , (1) 

where V is the vortex field operator whilst A, fi and ( denote constant param- 
eters. M is an integer bigger than N, to ensure that the energy is bounded 
from below. 

Vortex condensation corresponds to the spontaneous breaking of the mag- 
netic Z^ symmetry in In this case, domain wall defects can be formed, 
which for planar systems are, of course, one- dimensional objects. In order to 
be observable, a pair of heavy quarks must be interpolated by this kind of 
one-dimensional 'string', what results in a configuration energy which grows 
linearly with the pair separation distance pQ. 

A concrete realization of this idea had been previously introduced for 
the case of the SU(2) Georgi-Glashow model [2]. In the London limit, the 
relevant long-distance physics is dominated by the (unbroken) U(l) subgroup 
and, in that regime, the theory behaves as compact QED3, in the sense that 
it admits instantons. Far from their center, the instanton gauge fields along 
the locally unbroken U(l) direction look like Dirac monopoles (where one of 
three coordinates is the Euclidean time), and their action is UV finite since 
they have a finite size core. By using a singular gauge transformation, the 
original (monopole-like) instanton configurations can be also gauged along a 
fixed U(l) subgroup (unbroken by the vacuum) to be seen as Dirac monopoles 
in the whole space (with a corresponding Dirac string). This defines an 
effective Abelian theory (see also |H|) where the above mentioned finite size 
core may be interpreted as corresponding to a minimum length (i.e., an UV 
cutoff), a feature one should have, since compact QED is usually defined on 
a spacetime lattice. 
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Thus, as explained in |2j, a consequence of the compactness of the under- 
lying SU (2) theory is the existence of (finite action) instanton configurations, 
which for pure compact QED 3 are properly described by a Coulomb gas. The 
dual description is then given by the (Euclidean) Lagrangian density, 

C-duai = —^d^d^ + £/i 3 cos . (2) 

A comparison with the N = 2 model in (JJJ shows that (j2J) represents a 
vortex model where the global Z 2 symmetry is spontaneously broken (vortices 
condense), with the natural identification V = V exp(i(p/2) (see also [I]). In 
terms of this dual model, Polyakov obtained an area law for the Wilson loop 
in pure compact QED3, a signal of the confinement of electric charges. 

Several authors have afterwards discussed the effects due to the inclusion 
of parity preserving dynamical matter in the fundamental representation, a 
nontrivial modification of the model, since the symmetry is derived under 
the assumption that matter, when present, is in the adjoint representation. 

In jHlin]; the effect of massless fermions was taken into account by keeping 
just the quadratic part of their contribution to the effective action, a proce- 
dure which is justified by invoking an IR fixed point argument. As it is well 
known, when parity is not broken, the quadratic part of the fermionic effec- 
tive action contributes a nonlocal Maxwell term S n i to the effective action, 

Snl = \tj d 3 X U-^=U , (3) 

where / M is the dual of the field strength tensor corresponding to the gauge 
field A^, namely, = e^xdyAx, and £ is a dimensionless constant which, 
for a single fermionic flavour is £ = jg. In particular, in jH], based on the 
instanton anti-instanton logarithmic interaction implied by (jSJ), arguments 
in favour of a deconfined phase associated with instanton suppression were 
presented. 

In these references, compactness of the U(l) model was simulated by 
introducing instanton gas configurations (whose action is decreed to be finite) 
into the anomalous compact QED% model based on (jHJ). As a result, a dual 
description displaying a global Z 2 symmetry was obtained, 

S = ^d^i-d^dy + ^cos^, (4) 

which corresponds to an anomalous version of (j2J). 

However, as we will see, the computation of the Wilson loop based on the 
effective model (JIJ) would give a destabilization of the area law, instead of 
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the perimeter law characterizing deconfinement. This can be traced back to 
the global Z 2 symmetry of the anomalous dual action (jlj) which implies that 
the relevant source when computing the Wilson loop is indeed concentrated 
on a surface. 

On the other hand, in it had been argued that in the SU(N) Georgi- 
Glashow model with dynamical matter, the global becomes local, im- 
plying that the Wilson operator changes the vacuum inside the associated 
Wilson loop to a gauge equivalent state. Then, as discussed in [Zj , the overlap 
between the vacuum and the new state is expected to give a boundary effect 
and, depending on the localization properties of the model, a perimeter law 
would instead be obtained. 

In jH], compact QED 3 with dynamical matter was studied by using a 
bosonized version of the fermion sector. In this way, the local Z 2 theory 
invoked in [2j emerges naturally, with the bosonizing field playing the 
role of the Z 2 gauge field. Moreover, in this reference, the source which 
is relevant to compute the average of the Wilson operator is shown to be 
concentrated on the loop and, for massless fermions, arguments in favour of 
the perimeter law were presented. 

The point here is that, even in the massive case, the local Z 2 is still 
present. Thus, when the fermions become heavy, a destabilization of the 
perimeter law is obtained, instead of the expected area law. 

The question suggests itself about how, by tuning the fermion mass, dif- 
ferent low energy phases arise. When the dynamical fermion mass is large, 
fermions are almost decoupled and the area law (confinement) correspond- 
ing to the quenched theory should emerge. On the other hand, for light 
fermions a perimeter law characterizing deconfinement is expected. In this 
regard, more recently, a U(l) 3D compact lattice gauge theory having non- 
local interactions simulating the effect of the fermionic effective action has 
been studied (see jH]). Depending on whether a localized (representing mas- 
sive fermions) or nonlocalized interaction (representing massless fermions) is 
considered, the fitting of the Wilson loop using the lattice data prefers the 
area or the perimeter law, respectively. 

Our main proposal here is to discuss how the fermion mass is reflected 
in the corresponding low energy effective dual description of the continuum 
theory. Once the proper effective description is established, we will deter- 
mine the corresponding Wilson loop behaviour, analyzing confinement vs 
deconfinement . 

The structure of this article is as follows: in section |21 we start by analyz- 
ing the main properties of the Wilson loop in low dimensional noncompact 
gauge theories with and without dynamical fermions. This analysis will turn 
out to be useful when discussing some compactification scenarios in the pres- 
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ence of massless fermions. In section we will compare the Wilson loop 
representation implied by the global and local Z 2 dual descriptions, respec- 
tively. 

In particular, we will show that the local Z 2 symmetry amounts to the 
complete erasure of instantons, that is, the theory becomes effectively non- 
compact. 

In section HJ we analyze the effect coming from quasi- zero modes present 
in the zero charge instanton sector in compact QED with massless fermions, 
interpreting the result in terms of a binding interaction between opposite 
instanton charges. In this manner, we discuss the physical conditions that 
render the dual model with local Z 2 symmetry reliable, showing that the 
inclusion of massless fermions in compact QED 3 drive the system into a 
deconfmed phase. 

In sectional we study possible deviations from the local Z 2 (noncompact) 
scenario introducing a novel description when instantons are not completely 
suppressed: the instanton dipole liquid. In this manner, we will be able to 
discuss the stability of the deconfmed phase. 

Finally, we present our conclusions, comparing our approach with previ- 
ous ones discussed in the literature. 

2 Wilson loop in noncompact models 

A quite useful object for the study of confinement in gauge theories is (W(C)), 
the vacuum average of the Wilson loop operator W(C) corresponding to a 
(smooth) closed curve C, which (in d spacetime dimensions) is defined by 

W(C) = e ^c dx ^ . (5) 

It can also be written in an equivalent form by introducing the dual of the 
F^y field. For example, in d = 2: 

W(C) = e^ d2xxsf , (6) 

where f(x) = \^n V F^ v {x) and Xs(^) is the characteristic function of the plane 
surface E, enclosed by C. In d = 3 we have instead: 

W(C) = jf**'"'" , s„ = g5xn,, (7) 

where 5^, is a Dirac distribution with support on a spacetime surface £ = £(C) 
whose boundary is the curve C, and n M is the field of unit vectors normal to 
E, defined by the surface element dS^ = dS. 
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The vacuum average (...) is evaluated differently depending on whether 
the theory includes or not matter fields. In the latter ('quenched') case 

J DA ... e~ SG[A] 
= JVA e~ s o[A] ' ( 8) 

where Sg[A] is the gauge field action (we omit writing gauge fixing objects 
explicitly, since we assume them to be part of the measure VA). When 
charged matter is included, we have instead: 

JVA ... 

~ fVAe- s olA]-r F [A] ' ( 9 ) 

where IV [A] is the contribution of the matter (fermion) fields to the gauge 
field effective action: 

e -r F [A] = J v ^ v ^ e -s F m,A) ^ ( 10 ) 

with 

S F $,il>,A) = J d d xi)Vi) (11) 

and T> =<ft + % A" + m . 

In the two following subsections we evaluate the Wilson loop in non- 
compact QED2 and QED3 discussing the effect of introducing dynamical 
fermions. Let us begin with QED 2 , a warming-up exercise, where every- 
thing may be exactly calculated (at least when the fermions are massless). 

2.1 Noncompact QED2 

Let us first consider (W(C)) for the quenched noncompact theory: 

fVA e -S G [A)+ifsA 

W C » = SVAe-ScW ■ < 12 > 

where Sq = ^2 / d 2 xf 2 , and g is the coupling constant. We then introduce 
an auxiliary pseudoscalar field X(x) to write the equivalent expression: 

, , „ fVXVAe-^ x2+ ^f( x+ ^ , s 

JVXVA e -I a 2 x2+l If x 

Integrating A^ yields the constraint A + xs = K , where k is a constant. Its 
value is fixed to zero by using the boundary condition X(x) —> for x — > 00. 
Using this constraint into (|35|) yields 

(W(C)) = e -^f d2x ^ )2 = e -^^ d2x ^ . (14) 



6 



Finally, if Ay, denotes the area of S (the area enclosed by C) we have just 
shown that: 

(W(C)} = e~^ As , (15) 

as it should be, since we know that the model is confining. 

When a Dirac matter field is introduced, we have to add the effective 
action contribution, which for massless fields is given by: 

F F (A) = (27T)- 1 J d 2 xf^f. (16) 

An entirely analogous calculation to the previous one yields now the result: 

(W(C)) = e -^* EO " lxs , (17) 

where 

O = 1 - —d~ 2 . (18) 

7T 

For the particular case of a circular curve of radius R, the previous ex- 
pression may be exactly calculated. Indeed, by introducing the Fourier rep- 
resentation for the kernel of O^ 1 , 

0~ l (x v) = I ^ e iHx ~ y) - (19) 

V'V) J (27r) 2 e l + £> ^) 

we see that 

v 0-*v - I £k m)? (20) 

X(k) = J d 2 xxn{x)e^ = ^J^kR), (21) 
with Ji denoting a Bessel function of the first kind. Then 

, 2 



where 



In 



(W(C)) 



Ir* W?*L, m 

2 J k 2 + 9 -R 2 V } 



which can be exactly calculated: 
In 



(W(C)) = -TTg'R'hi^R)^^). (23) 

7T \/7l 



Here, I\ and K\ are modified Bessel functions. When R is very large (~/=-R >> 
1), we may use the asymptotic form : 

h{x)K x {x) ~ (24) 
— oc 



^^Next-to- leading terms are exponentially small. 
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to see that 



In 



(W(C)) 



7T2 



- — g X 2-kR 



(25) 



i.e., a perimeter law, a clear signal of the fact that the massless fermions drive 
the system into a deconfined phase. Of course, the area law is recovered if 
the g — > asymptotic form is taken for the Bessel functions. 



2.2 Noncompact QED 3 

As in the QED2 case, we first consider the quenched noncompact theory. 
Of course, here we have an expression for the Wilson-loop average which is 
formally equal to the one introduced in (|12j). except for the necessary changes 
due to the different number of spacetime dimensions. Namely, the gauge-field 
action Sq is now: 

while the term corresponding to the 'coupling' of the loop with A M , f cPxs^f^, 
deserves a closer look. To write it in a more explicit fashion, we first note 
that the 5% distribution may be given a more concrete form by introducing 
a parametrization of the surface bounded by the curve C: 

X^ : [0, 1] x [0, 1] - M 3 , (27) 

where X ^(a 1 , o~ 2 ) is a parametric representation of S. The two worldsheet 
coordinates a a , a = 1,2 are to some extent arbitrary, as well as their domain 
of definition. Nevertheless, 5^,, must of course be reparametrization invariant. 
Then we see that <5s is: 

5 E (x) = [ da 1 [ da 2 >/gl^j 5 {3) [x - X(a)] , (28) 
Jo Jo 

where g(a) denotes the determinant of g a b, the induced metric on the surface: 

g(a) = det [g ab (a)} , g ab {a) = ■ (29) 

The 'g(c) factor is of course required in order to have an invariant volume on 
the worldsheet, and as a consequence a reparametrization-invariant definition 
for <5s- To complete the construction of s M , we have to give an expression for 
hn(a), the normal to £ at the point parametrized by a. The surface element 
dS^(cr) is: 

ds,(a) = w§^§^f , (30) 



s 



which points in the normal direction, but is not normalized. Hence: 



\dS(a)\ 

and an elementary calculation shows that: 

\dS\ = a/ dS^dS^ = y/det(g t 

Thus: 



ui, ) 



(x) = da 1 da 2 VjR^^ - X(a)] h,(a) 



dX u 8X P 
dai d(J2 



'0 JO 

and using the explicit form of n„, we obtain 



(31) 



(32) 



(33) 



[ da 1 [ da 2 5 {3) [x - X{a)] e^ p 
Jo Jo 



(34) 



Now an auxiliary pseudovector field A^(x) is introduced in order to write an 
equivalent expression where only appears linearly in the exponent: 



(W(C)) 



VXD A e"V/ a m+*//m( a m+ s m) 



(35) 



(we have absorbed the normalization factor into the measure). Again, the 
integral over A^ yields a linear constraint: A M + s M = <9 M </>, where (ft is an 
(undetermined) pseudoscalar field. Then A can be completely integrated out 
by using that constraint, but there still remains an integration over 0: 



(W(C)) 



V(f) exp 



y / d 3 x (s, 



(36) 



It should be clear that the effect of the integration over is to erase any 
dependence of the average on the longitudinal part of s M . Indeed, 



In 



(W(C)) 



d 6 x s A1 (x)s A1 (x) 



+ -g 2 I d 3 x I d 3 y d ■ s(x)G(x - y) d ■ s(y) 



where 



G{x - y) 



d 3 k e iHx - y) 



1 



(2tt) 3 k r - 



Ati \x — y\ 



(37) 



(3J 
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The divergence of s M is obviously concentrated on the surface £, and there 
is a non-trivial (partial) cancellation between the two terms on the rhs of 
(|37|1 . To see that more clearly, we consider the simple (and relevant) case of 
a planar curve, contained in the plane £3 = 0: 

Sp(x) = 5^6(x 3 )xz(x) , (39) 

where x denotes the two coordinates on the x% plane. Since d ■ six) = 
5 / (x3)xe( x ), an elementary calculation shows that the contribution of the s 2 ^ 
term (which contains a <5(0)) is cancelled by a contribution from the second 
term, but there is a non-vanishing remainder: 



In 



(W(C)) 



g 



d 2 x 







AxE(x)[^G(x-y,i 3 )] I3=0 XE(y) • (40) 



Or, in terms of xs, the (two-dimensional) Fourier transform of xs, 



In 



(W(C)) 



d 2 k |k| 
(2vr) 2 T 



fe(k)| : 



(41) 



We note that the cancellation between two contributions in (jSZJ) gets rid of 
a divergent term proportional to the of s. 

To derive concrete results for (W), it is convenient to assume a particular 
form for the curve C (and its associated surface £). When C is a circle of 
radius R, using x(k) = J±(kR), we see that: 



In 



(W(C)) 



7T 



g 2 R 



AR 



du [Ji(u)f 



(42) 



where we have introduced an Euclidean UV cutoff, A. The integral over u 
can be exactly calculated; this yields 



In 



(W{C)) 



^g 2 R{KR)\F,(\\^\^-{KR) 2 ) 



(43) 



where 2-^3 is a generalized hypergeometric function. The last expression does 
not allow for a clean analysis of the dependence of the loop average with the 
distance, for long distances. This problem may be solved by using a milder 
form of regularization: let us consider the regularized function 



In 



(W(C)) 



1 



reg 



d 2 k Ikl 



(27T) 



IXE(k)P 



(44) 



with the obvious regularization parameter a. Note that, in order to keep 
the result dimensionless, we have had to introduce a constant fi, with the 
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dimensions of a mass. The regularization is removed by letting a — 1 — e, 
with e — > 0, and UV divergences emerge as poles in the complex variable e. 

Now the integral over k may be exactly performed for the analytically 
continued function, with the result: 



In 



(W(C)) 



7T 



r(|-f)r(i + f; 



J reg 



2 /i?( ^ Q 0Far(-f)r(2-f)- 



(45) 



Now we set a = 1 — e and collect the relevant terms when e — > 0, obtaining: 



In 



(W(C)) 



reg 



Ig* [±R + R hx(iMR)] (£-0). (46) 



The £ 1 term is of course divergent when the regularization is removed. A 
subtracted version of the average will then have the general form: 



In 



(W(C)) 



sub 



--g 2 [kR + R \n(nR)] 



(47) 



where the (finite) constant k is to be fixed by a renormalization condition. 
Of course, k may be absorbed into a different choice of the dimensionful 
constant //, and so we may use the expression 



In 



(W(C)) 



~ --9 2 [R Hfi'R 

- sub A 



(4* 



with a different \il ', to be fixed, for example, by setting the value of the average 
to zero at some small radius a, what yields the renormalized average: 



In 



(W(C)) 



(49) 



which is the result for the noncompact quenched case. Note that this result is 
consistent with the one obtained for a rectangular loop and using dimensional 
regularization ^U], where the static energy between point charges in d > 2 
is found to be: 



E(R) - E(a) = 

47T 2 

Taking the d — > 3 limit yields: 



,d-3 



r(^) (a 3 ~ d - R 3 ~ d ) 



E(R)-E(a) 



= 9 - In(^) , 

d^3 2n a 



to be compared with: 



E(R) - E(a) 



2 v a' 



(50) 



(51) 



(52) 
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for the circular loop, as follows from (|49|). 

Things are quite different when one includes fermionic matter into the 
game, since this adds the extra contribution already mentioned in equation 
© of the Introduction: T F (A) = S n i(A). 

Keeping only this contribution has the obvious advantage of having a 
quadratic gauge-field action. Thus, with just a few calculations we obtain 
now the Wilson loop average in terms of the Fourier transform of xs for this 
case also: 



In 



(W(C)) 



reg 



d 3 k \k\ 1+a 



(2 



7T 



k 2 



cat 

S k 



Ixu(k) 



(53) 



where we adopted the same analytic regularization as for the quenched case. 

Rather than evaluating directly this expression, we shall consider the 
results it yields in different limits. Of course, the £ — > limit is regular, in 
the sense that we recover the behaviour of the quenched case: 



In 



(W(C)) 



£^0 



(54) 



An interesting departure from this result arises when one considers the would- 
be 'large-£' limit, to be understood in the sense that only the IR-dominant 
contribution from the quadratic operator is kept. This is, indeed, justified 
for the calculation of the Wilson loop for a rectangular loop if one is only 
interested in the large- R case (in this case, g 2 R » 1. Then one has: 



In 



(W(C)) 



.1-a 



reg 



d 3 k \k\ a+1 



(2nf 



XE(k)| ; 



(55) 



which is of course independent of g. 

This integral can be directly evaluated, or its value can be deduced from 
known results corresponding to the d = 4 case, since the Feynman gauge 
propagator in this limit is identical to the one for d — 4, except from a 
constant factor. Indeed, one has: 



while 



U1 



(x\ 



-d 2 ) 



7T Z \X 



7T Z \X 



y\ 



(d = 3) 



(d = 4) 



(56) 



(57) 



and one is only interested in their values on a plane curve, so having a different 
number of components is irrelevant. For a rectangular loop, we obtain the 
result: 



In 



(W(C)) 



9_ 

7T£ 



r R 



(- - 1) Rg 2 » 1 



(5? 
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which shows a dramatic difference with the quenched case: the confining 
energy becomes a constant when massless fermions are included. 

Note that, as expected, (J49|) is associated with the logarithmic interaction 
between electric charges, typical of planar noncompact QED%, while the 
behaviour in ()58j) corresponds to a 1/R interaction between well separated 
electric charges. This comes about as, at large distances, (jo^j) is dominated 
by the action S n [(A) given in Q, and for planar systems, this kind of nonlocal 
action is known to imply a 1/R (Coulomb) static potential between electric 
charges (see [TT]). 

3 Compactification scenarios with fundamen- 
tal matter 

3.1 Instanton gas and QED%: Global Z2 symmetry 

Let us now consider the compact theory with a general nonlocal action. We 
shall consider the functional average: 

(W(C)) = J DiPAe-^ /(/+j)K(/+i)+i/s(/+i) , (59) 

where, to account for the neutral instanton background, we have made the 
substitution: 

Ux)^U(x)+Ux), (60) 

with 

jM = g m dy^S^ix-y), gg m = 2ir, (61) 

J-y 

where 7 is a set of Dirac strings associated with instanton anti-instanton 

(monopole-like) singularities located at the string endpoints. We have used 
a compact notation in the exponent, so that, for example, 

JfJCf = J a?xSyf^x)lC{x,y)f,{y) , (62) 

etcetera. Here, JC(x, y) is the kernel of an operator which, in the coordinate 
representation, is in general nonlocal. 

After introducing the auxiliary field A and integrating over the field, 
we obtain 

(W(C)} = [vjV<j>e-S KWHW+i/rt* . ( 63 ) 
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Now, if an integration over a dilute monopole plasma is supposed, 

N q a =±g m a 

= Y t N / N] { / dx (e i9mHx) + e - l9m ^ x) ) | , (64) 



N 

the following representation is obtained, 



(W(C)) = I V(f>e~ Id3:c ^ {d * +9mS)IC 1 C^+«m.)+2 fl g,Co«0] ^ (g5) 



where the field redefinition g m <f) —>■ <fi and Dirac's quantization condition 
g g m = 2n has been considered. 

In the local case, K.(x, y) oc S(x — y) and (jB3j) is just Polyakov's dual 
representation for the Wilson loop for pure compact QED%. According to 
Polyakov's analysis, the saddle point is given by a domain wall containing 
a 2tt discontinuity at S, coming from the source term. Therefore, when 
computing (<90 + g m s) 2 in (JBHj). the associated delta singularity in d<p cancels 
against the g m s = 2ir S^, n term. Because of the localization scale, away from 
E, the field is exponentially suppressed. Then, for a large surface, the domain 
wall leads to an almost uniform action density localized on S implying the 
area law for the Wilson loop, that is, linear confinement. Note that this 
analysis would also apply to the case where large mass fermionic matter is 
included, as in this case the fermionic effective action would contribute with 
just a local Maxwell term. 

On the other hand, when massless fermions are present, including the 
result from the one-loop contribution given by S n i (cf. ©) in the quadratic 
approximation, we see that 

K(x,y) = ( I |[l + -^]|y> 

At large distances, the second term in ()66|) dominates over the initial 
pure local Maxwell term. In this case, when disregarding the source terms 
in (jnnj), the anomalous sine-Gordon model considered in El (cf. (Jl}) is 
obtained. As in this model the relevant source to compute the Wilson loop 
continues to be concentrated on the surface £ (cf. (jHSJ)), or in other words, 
the Z 2 symmetry is still global, a perimeter law is not established in this 
framework, as required by a deconfining phase, but rather a destabilization 
of the area law due to the nonlocal kinetic operator. 
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3.2 Local Ztv and the Wilson loop 

In 0, by using the vortex condensation approach to the (2 + 1)D Georgi- 
Glashow model, the gauging of the symmetry in (JTJ) has been proposed, 
when matter in the fundamental representation is included. In this manner, 
a simple argument to discuss the Wilson loop behavior was presented. We 
review it here: In the absence of dynamical quarks, the (spontaneously bro- 
ken) Zfq symmetry is global. The VEV of W(C) is the overlap of the vacuum 
state |0 > with the state \S >, which is obtained from the vacuum by acting 
on it with the operator W(C). When acting on the vacuum, W(C) performs 
a Zn transformation at all points within the area E bounded by C. If the 
vacuum wavefunction depends on the configuration of non Z/v-invariant de- 
grees of freedom (the state in question is not Z^ invariant) the action of 
W(C) affects the state everywhere inside the loop. On the other hand, in 
the local theory with finite correlation length the overlap between the two 
states approximately factorizes into the product of the overlaps taken over 
the region of space of linear dimension of order of the correlation length I 

< 0|E >= n x < 0^ > , (67) 

where the label x is the coordinate of the point in the center of a given small 
region of space. For x outside E, the two states |0 X > and \S X > are identical 
and therefore their overlap is unity. However for x inside E the states are 
different and the overlap is therefore some number e -7 < 1. The number of 
such regions inside the area is obviously of order S/l 2 and thus 

< W(C) >= exp{- 7 ^} . (68) 

The VEV of W(C) then falls off as an area. 

On the other hand, when Zn is local, the action of a Wilson loop of 
finite size inside the contour corresponds to a gauge transformation. Thus in 
physical terms locally inside the contour the new state is the same as the old 
one and so the overlap between the two locally must be equal to one. The 
only nontrivial contributions to the overlap may come from the region close 
to the contour, thus (depending on the localization properties of the theory) 
a perimeter law is obtained in this case. 
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3.3 Bosonization and compact ificat ion of QED%: Local 
Z2 symmetry 

In jH], before compactifying QED 3 with dynamical fermions, the bosonized 
version of the matter sector was considered (see [T^l I13j). 

KfWi + ig J d^xA^J^ <-> Kb [A] + ig J d 3 x A tl e ilvp d v A p , (69) 

where Kp is the fermion action, and the U(l) fermion matter current is 
bosonized using a vector field A p whose gauge invariant action K B is given 
by the transverse Fourier transform of the fermionic effective action. The 
relation in (JBTJj) is exact, and refers to the equivalence of the partition func- 
tions associated with the fermionic field ip and the bosonizing gauge field A, 
respectively. As discussed in jH], the compactification process in this frame- 
work implies that the Z 2 symmetry present in (J2J) becomes local. In this case, 
the Wilson loop representation is (see jH]), 

lf~ - f d 3 x ( -VW+3 m s-27ri) 2 -^ 3 cos<£ }-K B [A] , 

W(s) = JVAV<j>e \^ yv9 ' V . (70) 

As the dual action is gauge invariant for any associated transforma- 

tion A — ► A — (l/27r)<90, the field in the correlator (|70|) decouples. In 
addition, parity symmetry implies that Kb [A] does not contain a Chern- 
Simons term; it only depends on the dual tensor edA. Then, transforming 
A — >• A + (g m /2ir)s, the exponent in (f7T)|) contains the term J d 3 x ^A 2 , plus 
the bosonized action which now depends on edA + (g m /2ir)eds. 

This led to a qualitatively different situation when compared with eq. 
(|65|). where the relevant source for the Wilson loop is concentrated on a 
surface. The local Z2 symmetry changes the relevant source to be e<9s| M = 
e fJ,l/p n p d I/ 5j:, which is tangent to the boundary of E and is concentrated there. 

This change in the dimensionality of the source indicates that, depending 
on the localization properties in (|70|l . a perimeter law could be observed in 
the phase where dynamical quarks are coupled. 

In order to evaluate the Wilson loop in the local Z 2 theory, gauging out 
the field 0, we obtain, 

= J_ J VAVAe~ f d3x lf 2+i I d?x{s-gA)edA~K B [A] ^ (y-g 
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and using the bosonization formula ([69)1 . this is precisely the Wilson loop in 
the noncompact model including the dynamical fermions, 



Then, we clearly see that the effect of the dynamical matter, via the 
associated local Z 2 symmetry, is that of gauging out the field <fi or, in an 
equivalent manner, completely erasing the effect of the instantons. That 
is, the theory becomes noncompact. In the next section, we will see that 
the complete erasure of instantons only occurs when the dynamical matter 
includes massless fermions, and certain conditions are satisfied, which settle 
the applicability of the local Z 2 scenario. In addition, the following remarks 
are in order: 

i) The resulting effect of the local Z 2 symmetry is similar to that ob- 
tained when introducing parity-breaking matter, although the physical 
reason is entirely different. In the parity-breaking case, the instantons 
are neatly erased because of topological reasons associated with the 
induced Chern-Simons term coming from the fermionic effective ac- 
tion (see [SI). Indeed, demanding gauge invariance in the presence of 
instantons, one discovers that there can be no instanton. 

ii) In the bosonization approach, the Z 2 local symmetry of the dual model 
is present even in the large mass limit. As a consequence, as discussed 
in [Bj, the relevant source for the Wilson loop is still concentrated on 
the boundary, and the expected area law is not re-obtained when the 
large mass limit is taken but a destabilization of the perimeter law. 

iii) Of course, the global Z 2 sine-Gordon model correctly describes the 
decoupling of large mass fermions and the corresponding confinement 
of the electric test charges associated with the Wilson loop. 

iv) On the other hand, according to our discussion in subsection when 
the effect of massless fermions is studied by means of the anomalous 
sine-Gordon model with global Z 2 , given in 13 E], a destabilization of 
the area law is obtained, instead of the perimeter law characterizing 
deconfinement . 




(72) 
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4 Zero modes and the induced instanton anti- 



instanton interaction 

From the discussion in section |3J the question suggests itself about how the 
fermion mass singles out the proper low energy description of the system. 

In this regard note that for large masses, /C ~ const., and the instanton 
interaction read from (}63|) (of course, disregarding the Wilson loop source s M ) 
becomes the usual Coulomb potential in three dimensions. Then, a dilute 
monopole plasma calculation is indeed justified, and this leads to the dual 
model represented by ()65|) . which displays a global Z 2 symmetry. Then, the 
relevant source to compute the Wilson loop is concentrated on the surface 
bounded by the loop; this together with the localization properties of the 
usual sine-Gordon model (/C ~ const.) lead to the area law. 

On the other hand, for massless fermions, as discussed in [Hj, if the trun- 
cation of the fermionic effective action is accepted, according to (J63|) with 

~ (pWz tlie instanton interaction would be logarithmic. This posses 
some reservations regarding the applicability of the dilute monopole plasma 
approach and the consideration of the anomalous Sine-Gordon model to com- 
pute the Wilson loop. This explains the destabilization of the area law dis- 
cussed in section 01 instead of the expected perimeter law, when considering 
this model. 

Then, to construct a reliable description of compact QED 3 with mass- 
less fermions, it is important to have a good understanding of the effective 
instanton interactions in this case. In particular, we should know how the 
energy 2 of an instanton anti-instanton pair depends on their distance. 

A concrete way to study that effective interaction between an instanton 
and an anti-instanton separated by a (spacetime) distance L, is to consider 
the function Q(L) 



where A L denotes the gauge field corresponding to the instanton anti-instanton 
pair, and its explicit form is of course dependent on the number of dimen- 
sions. Sf is the Dirac action for fields minimally coupled to A L . 

In what follows, we present the discussion of the two and three-dimensional 
cases. 

2 We follow the usual practice of calling energy to the quantity which in fact is an 
action. It may of course be regarded as an energy if one thinks of the Euclidean theory in 
d spacetime dimensions as a static theory in d + 1 dimensions. 




(73) 
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4.1 Two dimensions 



In two spacetime dimensions, we have for the Al gauge field the following 
condition: 

e^A L u {x) = 2irq[5 {2 \x - u) - 5 {2 \x - «)] (74) 

where q is the (integer) instanton charge, while u and v are the coordinate 
vectors corresponding to the instanton and anti-instanton locii, respectively. 

We then introduce bosonization, in order to calculate Q. Under this 
transformation, the fermionic current is mapped into the curl of a pseu- 
doscalar field 0: — > -^e^ u d u (j), while the fermionic action has a different 
form depending on whether the fermion mass is equal or different from zero. 
In the former case, we simply have a quadratic action for 0: 

Sf = \ j d 2 x (<90) 2 . (75) 

Thus, for the massless fermion case, 

Q(L) = -ln(e^ d2 ^)) = _ ^(e-'V^ff [*(«)-*(«)]) ? (76) 

where (...) denotes functional average in the bosonic theory, with the quadratic 
action (|75|). A straightforward calculation yields: 

Q(L) = 2q 2 \n\u-v\ , (77) 

where we have subtracted a contribution corresponding to the self-interactions 
(which is, of course, L-independent). In this regard, we note that the would 
be self-action of the vortex field due to the fermions, is divergent in the con- 
tinuum (of course, it will be finite if the vortices are given a finite core size): 



^-iVteqtK*)) _ oo ) m = 0. (78) 

In the massive fermion case, the bosonic theory is also exactly known, it 
becomes a massive sine-Gordon action: 

S F -> i J d 2 x[(d(f)) 2 + A cos(v / 4T0)] , (79) 

where A is a mass parameter. 
Of course, we still have: 

Q(L) = -ln( e -^ 9 ^-^) , (80) 
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where now (...) is evaluated with the action ([79)1 . This kind of correlation 
function is exactly known. Indeed, introducing the dimensionless variable: 
x = mL, and defining: 

w q (x) = -\L^-Q(L)] , (81) 

CLIj J x=mL 

it can be shown that w q may be expressed in terms of another function v q 

dyyv 2 g (y) (82) 
which satisfies the non-linear differential equation: 

< + -< = "7^4 «) 2 + «i - < + % ■ (83) 

y x h i _ v z h h x z i _ v z 

Rather than giving approximate numerical solutions for the equation 
above, we comment on the exact results that may be obtained in some par- 
ticular cases. We first note that, when the fermions are massive, the average 
for a single vortex is no longer zero \TE\ : 



e -*V4^(*)) = exp<j / - — ^ 7 - q z e- 2t } . (84) 



,m,„2 \ r°° dt 



2' l./ t 



sinh (qt) 



sinh 2 (t) 



This means that the would-be vortex self-action E = — \n(e~ l ^ i7Tq ^) , goes 
like E ~ — ln(am), where a denotes a small length scale, associated to the 
vortex size. This length scale defines the mass for which the self-action 
of the vortex is zero, which turns out to be m ~ 1/a, i.e., of the order of 
the cutoff. This is to be expected, since no quantum effect would come from 
those heavy states. 

Besides, one can also show [TH] that 

/ e -iV^ql(t>(u)-(j)(v)]\ _^ ^ e ~iVI^q4>^ e iV^q<f>^ ^ £ > oo . (85) 

In the massless case, we know that each factor is zero, and that the correlation 
function tends to zero as a power. In the massive case, however, we see that 
the vortex (q = 1) antivortex (q = —1) interaction energy coming from the 
fermionic sector tends to a finite constant when they are separated by an 
infinite distance: 

Q(L) 2 ln(am) , m ^ , L -»• oo . (86) 

We conclude that, in the massive case, the additional vortex interaction 
coming from the fermion sector is suppressed with respect to the logarithmic 
vortex interaction coming from the gauge field fluctuations. In other words, 
the effect of massive fermions is irrelevant. In the massless case, the induced 
vortex interaction is logarithmic thus renormalizing the effect of the gauge 
field fluctuations. 
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4.2 Three dimensions 



In the 2 + 1 dimensional case, we have already seen that the effect of large 
mass fermions is inducing a fermionic effective action of the (local) Maxwell 
form, simply renormalizing the instanton 1/L Coulomb interaction associ- 
ated with the gauge field fluctuations. In this case, when an instanton/anti- 
instanton pair is far apart, the fermionic effective action tends to a constant, 
or in other words, the fermion determinant does not suppress the instanton 
contributions. 

On the other hand, for massless fermions we do not have the exact so- 
lutions of the 1 + 1 dimensional case at our disposal. However, by using an 
approximated analytical treatment, which assumes that the monopole and 
antimonopole are far apart, one can find the effective interaction, even at fi- 
nite temperature [H| in the Georgi-Glashow model. This model is of course 
known to be a 'regularized' form of compact QED%. The T = result is 
found to be [T7j : 



where a is a matrix element obtained from the overlap between the zero 
modes corresponding to fermions in the background of a monopole and an 
antimonopole (with different centers). Explicitly, 



where /i is an IR cutoff. 

We note that in the m — > limit, at T = and for large L, the induced 
interaction coming from the massless fermion sector is not logarithmic, im- 
plying that the fermionic determinant is not suppressed as a power of 1/L. 
In other words, the fermionic effective action in (J3J) does not capture the in- 
stanton interaction at T = 0, as it would imply a logarithmic one. However, 
according to the discussion in we have the following situation: 

i) At high temperatures, dimensional reduction occurs. In the pure com- 
pact theory, the instanton flux lines are spread essentially on a 2D re- 
gion and this gives logarithmic interactions between instantons, which 
can be associated with a BKT transition. Above the critical tem- 
perature T c the instantons become suppressed forming dipoles. This 
corresponds to deconfinement at high temperatures. 

ii) When massless fermions are included, the critical temperature T c , above 
which the mean instanton distance is finite, is drastically lowered so 
that the dipole phase survives up to very low temperatures. 





a 



4n ln(/iL) , 
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iii) In this phase, for large distances, the fermionic determinant is sup- 
pressed as an inverse power of L: l/L 2N f. This means that above T c 
the anomalous quadratic contribution in (jHJ) captures the logarithmic 
behavior of the instanton interactions. As T c is exponentially sup- 
pressed, this means that the anomalous model can be used up to very 
low temperatures. 

iv) In the London limit, where the mass of the W± modes tends to infinity, 
the mean instanton distance in the molecules collapses to zero. Then, 
in the London limit, which corresponds to compact QED3, the theory 
including massless fermions becomes essentially noncompact. 

Because of i)-iv), the phase associated with compact QED 3 when including 
massless fermions, can be studied in terms of the noncompact version of the 
model. Taking into account the anomalous induced action in (JHJ), and accord- 
ing to the detailed calculation in subsection E21 this means a 1/R interaction 
between static charges and a Wilson loop perimeter law characterizing a de- 
confining phase. 

Note also that, as discussed in subsection I3.3| the noncompact model can 
be also rephrased in terms of a local Z2 symmetry. 

5 Instanton dipole liquid model: stability of 
the perimeter law 

The considerations above justify an effective dual model to describe the in- 
stanton phase above T c , taking as starting point the noncompact theory in- 
cluding the anomalous quadratic action in (J3J). This anomalous noncompact 
model represents the situation where the instantons are completely erased. 
Thus representing the Nf > 2 London limit of the Georgi-Glashow model. 
Deviations from this regime, as a large but finite mass for the W± modes 
or a small but nonzero fermion mass, can be analyzed by incorporating the 
formation of dipole configurations. 

We have already seen that the summation over instantons cannot be 
implemented as an anomalous instanton gas. We propose instead the sum- 
mation over configurations of the instanton dipole liquid type. That is, 



where M M is the magnetization field associated with the dipole instanton 
liquid and m = edM is the corresponding (topologically conserved) instan- 
ton current. We have included a minimal coupling between that conserved 
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current and the gauge field (the % is required in Euclidean spacetime), and 
S(M) denotes the action for the magnetization field, about which we shall 
say something below. 

The kernel O contains both the initial (pure) Maxwell term plus the 
nonlocal induced term in (JHJ), 

= l + (-^=, (90) 

note that, for large distances, the second term dominates over the local 
Maxwell term. 

For the sake of simplicity, we will assume that S(M) is a quadratic func- 
tional of the magnetization, 

S(M) = X - J d'xM^M^ (91) 

where Q is the kernel that defines the action for the magnetization field M. 
By symmetry reasons, it is reasonable to assume it to be a scalar function of 
the Laplacian. 

Of course, the integral over M and A can be performed exactly with the 
outcome: 

W(s) = e-l/^^o + n- 1 ) \± ; (Q2) 

where s 1 - is the transverse part of the current localized on the surface of the 
loop. 

Now, let us first assume the simplest (local) form for Q: Q = 1/a 2 , where 
a is a real constant. This represents a system where large magnetizations 
are suppressed, so that it describes a basic feature of the phase where the 
instanton dipole liquid is formed. In particular, note that in the a — > limit, 
the magnetization is completely eliminated and the noncompact situation is 
re-obtained. In the general case, when a finite a is considered, we see that: 



ll+a 



L -I reg 2 J (27r) d k 2 1 + £S_ _|_ fl 2 

(93) 

where we adopted the same regularization as in the noncompact case. 

This may be easily rearranged to look identical to the analog noncompact 
expression: 



In 



|l+a 



i 1 f rl 3 k Ikl 

W^L s -2^° J Wf Vl7^ fe(k)|2 ' (94) 
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with g a = g/y/1 + a 2 . Therefore, in this case, there is no departure from the 
perimeter law obtained in the noncompact case, except for a finite renormal- 
ization of the coupling constant. 

This comes about since the object that determines the behaviour of the 
Wilson loop average is (O + Q" 1 ) , so it is clear that the IR (long distance) 
behaviour of the loop might only be changed, with respect to the one dictated 
by O, if Q" 1 produced a stronger divergence for small momenta. Otherwise, 
there are not relevant changes. 

For this reason, if the above mentioned model is improved by including 
correlations in the magnetization at different points, Q — 4y (— d 2 ) + 1/a 2 , 
as the infrared behavior of O continues to dominate over the effect of the 
magnetization, this model will also display a perimeter law. 

An example where the IR behaviour may be changed corresponds to the 
a — > oo limit in the above expression for Q, that is, fl = \ (— d 2 ). Consider- 
ing the usual circular region in the x 3 plane, this leads, in Fourier space, to 
the following expression for the Wilson loop average: 



In 



(W(C)) 



d 3 k 



Xs(k) 



1 



+ 



(95) 



fc 2 



where xs is the (two dimensional) Fourier transform of the characteristic 
function for the circle. Using then the explicit form for the kernels O and Q, 
we find: 



In 



(W(C)) 



±R 2 
2 



dkn 



d\k\ |k| 



J 2 (\k\R) 



k 2 + E,g 2 k + k 2 



(96) 



which is UV divergent. Including the same kind of analytic regularization 
we applied in the free case, we see that: 



In 



reg 



(W(C)} 
or, by a change of variables 



1 /*+oo POO 

— R 2 (/j,) l - a / dk / d|k| |k| 

2 J-oo Jo 



J 2 (\k\R) 



k 2 + £g 2 \k\ + K 2 ' 
(97) 



In 



(W(C)) 



reg 



X 



--R(fiR) 



l-Q 



X 



111 



dL 



I 2 + ig 2 \l\R + k 2 R 2 



d\l\ 



(98) 



It should be evident (from the discussion in the free case) that the integral 
over shall produce a pole in e = 1 — a, and hence the result, for long 
distances, will be essentially RlnR. 
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Of course, the destabilization of the perimeter law is expected in this case, 
as an action for the magnetization based on the operator Q — \ {—d 2 ) would 
also represent a destabilization of the instanton dipole liquid: for homoge- 
neous configurations there would be no cost to increase the magnetization, 
or equivalently, the instanton anti-instanton distance. 

Conclusions 

In this article we have studied the mechanism by which the mass of (parity- 
preserving) fermions singles out the proper low-energy dual description of 
compact QED 3 , explaining whether the resulting phases are confining or 
deconfining. To that end, we have discussed the instanton anti-instanton 
interactions. 

For a large fermion mass, the quadratic approximation to the fermionic 
effective action was found to be reliable, giving a local Maxwell term which 
renormalizes the pure-gauge action. In this case, the instanton anti-instanton 
interaction obeys the usual Coulomb behavior. Then, a dilute monopole gas 
approximation can be considered to implement the compactification of the 
model, and the usual Polyakov dual model (with renormalized parameters) 
is obtained. This model possesses a global Z2 symmetry, so that the relevant 
source to compute the Wilson loop is concentrated on the surface bounded 
by the loop. Therefore, the area law characterizing confinement is obtained 
in this case. 

On the other hand, it is known that when two or more flavours of mass- 
less fermions are considered 3 , above an exponentially suppressed critical 
temperature T c , the instantons bind to anti-instantons to form dipoles. 

In the Georgi-Glashow regularized version of compact QED 3 , the aver- 
age instanton anti-instanton distance in the molecule is finite; however, in 
the London limit where compact QED3 is approached, the average distance 
collapses to zero. Then, in this situation, the effect of the massless fermions 
is that of completely erasing the instanton configurations, thus leading to 
an effectively noncompact gauge model. This model can be characterized by 
the particular suppression of the fermionic determinant, due to quasi zero 
modes, when the instanton and anti-instanton are separated by a large dis- 
tance L. As this suppression goes like an inverse power of L, it amounts 
to an attractive logarithmic instanton interaction, which is captured by the 
behavior arising from the quadratic truncation of the fermionic determinant. 

Note that changing the number of flavours requires, in our notation, to change the 
value of the constant £. 
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Summarizing, in the above mentioned conditions, the presence of massless 
fermions in compact QED% render the theory essentially noncompact, and 
hence the anomalous effective action induced by the massless fermions is a 
good starting point to study the corresponding phase. We have presented a 
careful calculation of the Wilson loop in this model to show that the perimeter 
law is indeed satisfied, and the corresponding phase is deconfining. This 
outcome is similar to that observed in compact QED 3 with parity breaking 
matter where the instantons were also erased, but due to topological reasons 
associated with the induced Chern-Simons term. 

On the other hand, for a large mass of the W± modes in the Georgi- 
Glashow model, but outside the London limit, or for nonzero but very small 
fermion masses, the average dipole size is expected to be finite. For this rea- 
son, we have also studied deviations from the noncompact situation by con- 
sidering a model for the instanton dipole liquid formed in this case, showing 
that the perimeter law is preserved for a reasonable class of actions charac- 
terizing the magnetization field for a liquid of stable instanton dipoles. 

Finally, with regard to the symmetry properties characterizing the low 
energy dual description of QED 3 , we have also shown, via bosonization, 
that the erasing of instantons can be associated to a situation where the 
initial global Z 2 symmetry in the pure gauge theory becomes local when 
massless fermions are included. This is to be compared with an anomalous 
dual model for compact QED 3 with massless fermions, previously discussed 
in the literature, possessing a global Z 2 symmetry. Here, we have shown that 
such a global Z 2 scenario would give a destabilization of the area law, instead 
of the expected perimeter law associated with deconfinement. 

Then, by tuning the fermion mass, from very large to very small values, 
there is a confinement /deconfinement transition which cannot be described 
by means of a single dual low energy effective model. The associated low 
energy descriptions are essentially different, displaying a global or local Z 2 
symmetry, respectively. 
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